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Abstract 

We study the dimensions of non-chiral operators in the Veneziano limit of M = 1 
supersymmetric QCD in the conformal window. We show that when acting on gauge- 
invariant operators built out of scalars, the 1-loop dilatation operator is equivalent to 
the spin chain Hamiltonian of the ID Ising model in a transverse magnetic field, which 
is a nontrivial integrable system that is exactly solvable at finite length. Solutions 
with periodic boundary conditions give the anomalous dimensions of flavor-singlet 
operators and solutions with fixed boundary conditions give the anomalous dimensions 
of operators whose ends contain open flavor indices. 



1 Introduction 



Conformal field theories in four dimensions are both ubiquitous and have many potential 
phenomenological applications (e.g., [1-29]). While much progress has been made towards 
unraveling the structure of theories with lots of supersymmetry such as M = 4 SYM, 
strongly-coupled CFTs with M = 1 or no supersymmetry remain fairly elusive. In particular, 
phenomenological applications often depend crucially on the dimensions of unprotected 
operators, and very few tools are generally available to learn about their dimensions. While 
some progress has been made recently towards placing general bounds on unprotected 
quantities [30-35], the situation is still unsatisfactory and it is important to develop new 
techniques to learn about 4D CFTs with little or no supersymmetry. 

In order to gain a better understanding of these theories it is very useful to find limits in 
which simplifying structures emerge. One such limit is that of weak coupling A 1 (possibly 
in a dual frame), where perturbation theory is valid. Another is a large- N limit [36], where 
the theory becomes planar, and operator dimensions and correlation functions factorize 
according to large- N counting. In the present work, we will focus on M = 1 supersymmetric 
QCD with gauge group SU(Af c ) in the conformal window |iV c < Nf < 3N C [37]. This theory 
has a large- N limit in the sense of Veneziano [38], taking N c — > oo while holding Nf/N c 
fixed. In the Veneziano limit it is additionally possible to take Nf/N c arbitrarily close to 
either edge of the conformal window, where either the electric or magnetic description of the 
theory has a perturbative Banks-Zaks fixed point [39] and is weakly coupled. 

The remarkable lesson of the AdS/CFT correspondence [40-42] is that large- N CFTs 
at very strong coupling A ^> 1 can be described in terms of perturbative string theory in 
AdS backgrounds. Moreover, the low-dimension spectrum of such theories can often be 
accurately described by effective field theories in AdS. Unfortunately, this limit is typically 
only accessible in theories with an exactly marginal coupling, such as = 4 SYM. By 
contrast Af = 1 SQCD flows to an isolated fixed point, where the 't Hooft coupling is 
generically expected to be 0(1). A quantitative string-theoretic description of these fixed 
points is currently lacking, and any such description is generally not expected to admit a 
traditional effective field theory limit [43, 44]. 

The last decade has also seen remarkable progress in understanding certain special 
theories (including M = 4 SYM and some deformations) at arbitrary values of the 't Hooft 
coupling. This is due to the property of integrability, which roughly means that the theory 
has an infinite number of hidden symmetries not manifest in the Lagrangian description 
(see [45] for a review and many references). Some of the first hints for integrability in 
M = 4 SYM came from calculations of the action of the 1-loop dilatation operator on closed 
sectors of operators built out of scalar fields. In [46] it was shown that the 1-loop dilatation 
operator acting on the SO (6) subsector of scalar operators (closed at 1-loop) is identical to 
the Hamiltonian for an integrable SO (6) Heisenberg spin chain, which may be diagonalized 
using a Bethe ansatz. This investigation paved the way to discovering integrability of the 
full 1-loop dilatation operator of the theory [47-50], and was soon generalized to various 
sectors of the gauge theory at the higher-loop level (e.g., [51-62]). 
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While it would be quite surprising to find higher-loop integrability in a wide range of 
theories, the success of the above program suggests that potentially interesting structures 
might lie hidden in other 4D CFTs. Moreover, it is not yet clear precisely what role 
supersymmetry plays in creating such structures. We believe that it is important to gather 
more data about how and when simplifying structures can emerge in theories qualitatively 
different from A^ = 4 SYM, particularly those with potential phenomenological relevance. 

Thus, in the present work we undertake a similar investigation to that of [46] for the 
Veneziano limit of M = 1 SQCD. We will show that when acting on gauge-invariant 
operators built out of scalars, the 1-loop dilatation operator (in the electric Banks-Zaks 
limit, Nf ~ 3iV c ) is identical to the Hamiltonian for the 1-dimensional Ising model in a 
transverse magnetic field. Here each site on the spin chain is occupied by a flavor-contracted 
(gauge adjoint) "dimer", either QQ^ or Q^Q. This is perhaps one of the simplest examples 
of a nontrivial integrable model, and can be mapped to a system of quasi-free fermionic 
excitations through a combination of Jordan- Wigner and Bogoliubov transformations. We 
exhibit the exact solutions at arbitrary finite length for both periodic and fixed boundary 
conditions, which give the 1-loop anomalous dimensions of operators containing both closed 
and open chains of scalar quarks. 

An important difference between this sector of operators and the sectors typically studied 
in M = 4 SYM is that the ground state is non-BPS, and receives corrections at each 
order in perturbation theory. Indeed, simply identifying the 1-loop ground state in terms 
of elementary fields is a nontrivial task. By contrast, in M = 4 SYM one can consider 
sectors of operators where the ground state is BPS and its form and dimension is protected 
to all orders. In such cases, the superalgebra preserving the ground state acts linearly on 
excitations above it [63], which is helpful in solving the dynamics and making all- loops 
predictions for the dimensions of operators [58, 59, 64-67]. We have not found a similar 
picture in A^ = 1 SQCD, and it is not yet clear how one would make precise conjectures 
about the spectrum deeper into the conformal window even if integrability were to persist. 

Another notable difference is that excitations above the ground state in M = 1 SQCD 
are highly nonlocal collective modes, in contrast to the local magnon excitations typically 
studied in M = 4 SYM. We expect that similar nonlocal excitations might be present in 
other theories with a Veneziano limit, such as Af = 2 SQCD with gauge group SU(iV c ) 
and Nf = 2N C flavors. The 1-loop dilatation operator of this theory was recently explored 
in [68, 69], with intriguing but inconclusive results about its integrability properties. We 
hope that the present study might be useful in better understanding the structures that 
appear in this and related theories. 

This paper is organized as follows. In section 2 we discuss the Veneziano limit of SU(Af c ) 
SQCD in the conformal window. In section 3 we derive the form of the 1-loop dilatation 
operator acting on operators built out of scalar fields, and show that it is equivalent to 
the Hamiltonian of the ID transverse field Ising model. In section 4 we review the exact 
solution to this system at arbitrary length for both periodic and fixed boundary conditions, 
and discuss the resulting spectrum of operator dimensions. We conclude in section 5. 
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2 SQCD in the Veneziano Limit 



We will be primarily interested in M = 1 SQCD with gauge group SU(iV c ) and Nf flavors 
of quarks, Q ai and Q la (i,i = 1...JV/, a = l...N c ). In the range |JV C < N f < 3N C 
the theory is believed to flow to an interacting conformal fixed point [37]. The theory has 
anomaly-free global symmetries SU(Nf) L x SU(Nf) R x[/(l) B x U(1)r, where Q transforms 
as (N f , 1, 1, 1 - N c /N f ) and Q transforms as (l,N f , -1, 1 - N c /N f ). 

The Veneziano limit [38] corresponds to taking N c — > oo while holding Nf/N c fixed. In 
this limit the theory contains a single continuous parameter e = 3N C /Nf — 1, and when 
e < 1, the electric description is weakly coupled and one can reliably perform calculations 
in perturbation theory [39]. The fixed-point value of the 't Hooft coupling then has the 
expansion [70] 1 

q 2 N c 1 o 9, s , . . 

One can also compute the dimensions of gauge-invariant operators order by order in e. A 
simple example is the chiral meson operator (QQ)], whose exact dimension is controlled by 
the super conformal U (1)r symmetry. In this case we have = \Rqq = 2 — e. Note that 
the baryon operators decouple in the large N c limit, and will play no role in our discussion. 

An important simplification that occurs in the N c , Nf — > oo limit is that computations 
in perturbation theory are dominated by planar diagrams. Since N c and Nf are on the same 
footing when it comes to large-A^ counting, one must keep track of both color lines and 
flavor lines in order to decide if a diagram is planar. As a consequence, the usual large- N 
factorization story is modified. In the Veneziano limit, operators factorize into products 
of "generalized single-traces," namely strings of fields that have adjacent color and flavor 
indices contracted (up to a possible overall trace). Focusing on operators built only out 
of the scalars Q and Q for simplicity, the basic building blocks of generalized single-trace 
operators are the flavor-contracted color adjoints X = QQ^ and Y = Q^Q. From these, one 
can construct the flavor-singlet gauge-invariant operators 

Tr(XY...X) = Q ai Q j!ib Ql d Q Jc ...Q dk Q j!ka , (2.2) 

as well as the flavor adjoint and bi-fundamental gauge-invariants 

(Q*XY...XQ) i j , (Q^XY...XQ% (QXY...XQ)], (QXY...XQ% (2.3) 

Other operators built only out of scalars can in general be written as a product of the 
operators in Eqs. (2.2) and (2.3), and will be referred to as "generalized multi-trace" opera- 
tors. They have dimensions that are equal to the sum of the dimensions of the constituent 
single-traces. 

Note that operators that are constructed by contracting left flavor indices with right flavor 
indices, such as a chain of chiral mesons Q ia Q a jQ^ b Qbk ■ ■ ■ Q lc Qd, are generalized multi-trace 

1 Here we give the expansion for the coupling as defined in the DRED renormalization scheme [71]. 
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operators. Their correlation functions and dimensions factorize into those of their constituent 
generalized single-trace parts. 2 If we further insert an "impurity" into the chain, for example 
by replacing one of the mesons with Q^Q, then the only nontrivial mixing in the anomalous 
dimension matrix will occur within the generalized single-trace operator that contains the 
impurity, QXQ. 3 Consequently, it is unclear how to realize in a nontrivial way the M = 4 
picture of a BPS ground state sprinkled with impurities, at least using only scalar quarks. 

The operators in Eqs. (2.2) and (2.3) have classical dimensions that are equal to the 
number of constituent quark fields. At the quantum level, operators carrying the same 
global symmetry charges can in general mix with each other, and the eigenvectors of this 
mixing acquire anomalous dimensions that are functions of e. At the 1-loop level (leading 
order in A ~ e) the above operators mix only with each other and form a closed subsector. 
At higher loops they can also mix with operators containing fermions, gauge fields, and 
derivatives. In the following sections we will focus on understanding the dimensions of these 
operators at the 1-loop level. 



3 The 1-Loop Dilatation Operator 
3.1 Feynman Graphs 

Let us first briefly recall how one computes the dilatation operator in perturbation theory. 
Consider a collection of bare operators Of 1 with two-point functions that are canonically 

normalized at zero coupling (Of ) \x)O^°\0)) a=o = 5ijX~ 2A * ■ For us, these will be the 
generalized single-trace operators of Eqs. (2.2) and (2.3). At nonzero coupling, their two- 
point functions have the form 

x^\of\x)Of\Q)) = 5 ij -\V ij (± + \og(v 2 x 2 ) + ..)j+0(\ 2 ), (3.1) 

where we have regulated divergences with dimensional regularization. 4 Now, renormalized 
operators Oi with well-defined scaling dimensions A, are linear combinations of bare opera- 
tors Oi = Z^Oj. The Aj have a perturbative expansion of the form Aj = Af ) ' ) -|-A7j-|-0(A 2 ), 
so that two-point functions of renormalized operators can be written 

a^Vi(*)0}(O)> = dyaT 2 ^-^) 

= 5ij - \%6ij log(x 2 ) + 0(A 2 ) 

= (ZZ% - X(ZVZ% Q + \og(n 2 x 2 )) + 0(A 2 ), (3.2) 

2 This can be verified directly by drawing Feynman diagrams with double-line notation as in [68]. With 
only the interactions of J\f = 1 SQCD, it is impossible to "fill in" a region of a diagram between left and 
right flavor lines. 

3 In other words, the dispersion relation for such an impurity will be independent of momentum. 
4 The reader should not confuse the dimensional regularization parameter e with the expansion parameter 
e = 3N C /Nf — 1 defined in the previous section. 
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Q Q g fia Qai g fia Qaj Q m Q ]lh Qai Q j 




Figure 1: 1-loop Feynman graphs contributing to logarithmic divergences in two-point functions of 
generalized single-trace operators in the planar limit. Dotted lines represent the auxiliary field D 
in the vector multiplet. All other relevant 1-loop graphs are obtained by replacing Q,Q^ Q\ Q 
on one or more scalar lines. 



where in the last line we have substituted Eq. (3.1). Matching the 0(A°) term and 1/e pole 



above, we find Z = U — ^UV + . . . , where U is a unitary matrix, and ". . . " represents 
scheme-dependent finite terms and general terms of higher order in A. Finally, matching 
the coefficients of Xlog(x 2 ) gives 7 = UVW, so that the 1-loop anomalous dimensions are 
simply the eigenvalues of V. 

Returning to SQCD, we have the task of computing 1-loop corrections to two-point 
functions of bare operators Of 3 ^ and extracting V - - the coefficient of their 1/e poles. 
In the planar limit, these corrections come from two sources. The first is wavefunction 
renormalization for individual fields making up Of\ for example graphs 1(a) and 1(b). The 
sum of these graphs turns out to vanish in Feynman gauge (where the Euclidean gluon 
propagator is simply g^/k 2 ). Hence we choose to use Feynman gauge in what follows. 

The second contribution to V is from interactions between pairs of neighboring fields 
111 Of\ Such pairs can either be gauge-contracted (e.g. graphs 1(c) and 1(d)) or flavor- 
contracted (e.g. graphs 1(e) and 1(f)). For gauge-contracted pairs, we must normalize each 
graph by dividing by the tree- level two-point function N c A 2 xy , where 

r(i - e) 1_ 

47T 2 ~ e \x — y\ 

is the scalar propagator in d = 4 — 2e dimensions. We find, 

1(c): -j^9%^{T A T A )^J d d z^ x A% = -5^5^+ finite, (3.4) 

1(d) : j^g%5]TT{T A T A )^J d d z d d w {A xz %A zy )( gflv A zw ){A x J»A wy ) 

= 5$^+ finite, (3.5) 

where A = is the 't Hooft coupling. 

The divergent parts of Eqs. (3.4) and (3.5) exactly cancel, so that nearest-neighbor 
interactions between gauge-contracted Q^Q pairs vanish at 1-loop. However, if we replace 
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one of the Q's with Q\ then the D-term contribution 1(c) flips sign, while 1(d) is unchanged. 
Consequently, nearest-neighbor interactions between gauge-contracted pairs take the follow- 
ing form at 1-loop: 



v; 



augc contr's 



/ 


gtg gtgt qq qq] \ 







QQ 


-i 


Q ] Q ] 


-i 


\ Q ] Q 


/ 



x AI 



flavor • 



(3.6) 



For interactions between flavor-contracted pairs, we must divide by the tree-level two- 
point function NfA 2 xy . The graph 1(f) vanishes by antisymmetry of the gauge-scalar-scalar 
coupling. The remaining contribution 1(e) is identical to 1(c), up to overall factors. Once 
again, the sign flips if we replace Q <-> Q\ so we find the following contributions to V from 
flavor-contracted neighbors: 



^flavor 



contr's 





gg f g f g 




i -i 




-i i 



N 



x 



A x 2(T A ) b a (T 



2N, 



(3.7) 



Whenever our flavor-contracted pair is part of a chain of more than two fields, the factor 



■5 a 5^ is equivalent to the identity (I 



\cb 
;auge ) a d 



5^5 d in the planar 



limit. However, there is an important exception for two-field operators, where each gauge 
generator is traced-over, giving Tt(T a )Tt(T a ) = 0. In this case, the naively subleading piece 
— jf5 b a 5d contributes at the same order in N c . 

From here, it's easy to read off the dimensions of operators containing two fields. Con- 
tributions from flavor contractions vanish, so that the 1-loop anomalous dimension matrix is 
given simply by Eq. (3.6). The mesons QQ have 1-loop anomalous dimension —A, consistent 
with their superconformal i?-charge. Meanwhile, both Q^Q and QQ^ have vanishing 1- 
loop anomalous dimensions, consistent with the fact that the baryon current Q^Q — QQ^ is 
protected to all loops, while the Konishi operator Q^Q + QQ^ has dimension equal to the 
slope of the /3-function at the conformal fixed point, which has its first nonzero contribution 
at two loops [72]. 



3.2 Spin Chain Hamiltonian for General Operators 

To discuss more general operators, it is useful to introduce the notion of operators as spin 
chain states, with V acting as a Hamiltonian. Generalized single-trace operators are built 
out of strings of gauge adjoints X = QQ} and Y = Q^Q, which we will call up states |t) and 
down states respectively. In this language, we may write 

^gauge contr's = ^ ® (J* - I ® I) , (3.8) 
^flavor contr's = 2JV"^ 
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where a x and a z are the usual 2x2 Pauli matrices. Then traces Ty(XYX . . . ) of length 
L correspond to states | tit • • •) in a circular spin chain with L sites. For such states, the 
Hamiltonian is simply a sum of Eqs. (3.8) and (3.9) for each site, 

L - 1 A N 

Closed = Y J 2 { <<^- l)jr T,W^ 1 -<) 

n=0 n=0 c 

N f -N A ^ / N f \ 

= ^^ 2 L ff A +1 -^ n , (3-10) 

c n=0 V c / 

where a z L = <Tq. Note that this Hamiltonian commutes with the parity operator P = rin=o a n-> 
and can be diagonalized separately on each of its eigenspaces P = ±1. 

Similarly, we can think of the flavor adjoint and bi-fundamental operators in Eq. (2.3) as 
states in an open spin chain, with nondynamical sites at each end. For example, Q^XY . . . XQ 
and QXY . . . XQ correspond to |y:T4 ••• T : T) an d ■■■ T : T); respectively, where we 

have separated off the nondynamical sites with colons. There is no flavor contraction at a 
nondynamical site, so the Hamiltonian for an open chain of length L + 2 is 

-f£«Y (3.1D 

c n=l / 

In the next section we will proceed to diagonalize Eqs. (3.10) and (3.11) in order to 
obtain the 1-loop spectrum of operator dimensions. Before we do so, let us mention that the 
generalization of these results to SQCD in the conformal window with an SO(iV c ) [37, 73] or 
Sp(iV c ) [74] gauge group (in the Veneziano limit) is entirely straightforward. In both cases, 
the only difference is that the states with odd parity are projected out (i.e., they vanish 
identically), but otherwise the Feynman diagrams and spectra of 1-loop operator dimensions 
are the same. 



4 Solving the Spin Chain 

4.1 Closed Chains 

Eq. (3.10) is the Hamiltonian for a 1-dimensional anti-ferromagnetic Ising model in a trans- 
verse magnetic field. Remarkably, this is a textbook example of an integrable model [75, 76], 
and we will be able to write down its energy levels and eigenstates exactly. The following 
solution is standard, but we include it here for completeness. 



n=0 n=l ' 

, N f -N A A /A 

c \n=0 
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We will study the Hamiltonian 



L-l 



n=0 

for a magnetic field h acting on a circular spin chain of length L. Since our Hilbert space is a 
tensor product of two possible states for each site, it's tempting to think of it as a fermionic 
Fock space, with creation and annihilation operators = h(Cn ± io~*) f° r each position 
n. This is not quite correct, since these operators commute rather than anticommute at 
different sites {a^,a^} ^ 0. The way around this problem is to perform a Jordan- Wigner 
transformation [77], defining fermionic creation and annihilation operators: 

n—l n—1 

-t = 



n<^n, C n EE H a - a -' ( 42 ) 



m=0 m=0 

which now satisfy canonical anticommutation relations 

{4' C m } 0~nrm {c n , C m } {4' 4i} 0. (4-3) 

These new creation and annihilation operators are nonlocal on the spin chain, and we can 
think of them as creating and destroying fermionic solitons. 

In these new variables, we have 

L-2 L-l 

H = ^(4 + cJ(c^ 1 - Cn+1 )-P(4_ 1 + CL-i)(4- C o)-/ i ^(2 C t Cn -l) (4.4) 

n=0 ra=0 

where P = Yln=o a n = ex P i}^ S n 4 c n) — is the parity operator. Within each parity 

eigenspace, the Hamiltonian becomes simply 

L-l L-l 

H = ^(4 + c n )(4 +1 -c n+1 )-/i^(24c n -l), (4.5) 

n=0 n=0 

where when P = — 1 we must impose periodic boundary conditions cl = cq, and when 
P = 1 we must impose antiperiodic boundary conditions cl = —Cq. 5 Our Hamiltonian is 
now translationally invariant and quadratic in creation and annihilation operators, so it can 
be diagonalized via a Fourier transform and Bogoliubov transformation: 

H = cos k + 2h)c\ck — i sin k [ c^_ k c\ + c_fcCfc j + L/t (4.6) 

k 

5 Note that states with P = 1 are exactly those created by an even number of c's or c*'s, so the states 
themselves will have periodic boundary conditions and can make sense on a circular chain. 
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where b k and bk are new canonically normalized creation and annihilation operators, and 

e(Jfe) = 2y/h 2 + 2hcas{k) + 1 (4.8) 

is the dispersion relation for the free fermionic quasiparticles created by b\. Because of the 
boundary conditions imposed by parity, the quasimomenta must take the values 



2rnir p i 

k = { { 2m L + l)n p Z ! ( 4 - 9 ) 



for m = 0, 1, . . . , L — 1. 



4.1.1 Examples: Closed Four-Field and Six-Field Operators 

Specializing to the case of interest, we can now write 

V closed = XL + X^ (b]jb k - l - J v / 10 + 6cos(A;), (4.10) 

where we have used h = jf ~ 3 near the weakly-coupled end of the conformal window. In 
addition to having quasimomenta obeying the quantization conditions in Eq. (4.9), states 
corresponding to traces Tr(XYX . . . ) must be invariant under shifts of the spin chain — in 
other words they must have vanishing total momentum modulo 2ir. 

As an example, let us consider four-field generalized single-trace operators, spanned by 
Tr(XX), Tr(yy), and Tr(XY). The two-site spin chain has two eigenstates of odd parity, 

bl\0) and 6j|0), (4.11) 

but only the latter has vanishing momentum. It corresponds to the operator Tr(XX) — 
Tr(yy), with anomalous dimension 

^2 + + 6cos(0) - ^10 + 6cos(tt)^ A = 3A. (4.12) 

Meanwhile, we have two eigenstates of even parity, both of which have vanishing momentum, 

|0) and &L&t|0). (4.13) 

2 2 

These have anomalous dimensions 

^2 ± + 6 cos(tt/2) ± + 6 cos(3tt/2)^ A = (2±v / l0)A, (4.14) 

and correspond to linear combinations of Tr(Xy) and Tr(XX) + Tr(yy). 

As another simple example, let us repeat the above analysis for six-field generalized single- 
trace operators, built out of linear combinations of Tr(XXX), Tr(XXy), Tr(Xyy), and 
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Tr(YYY). Using the L = 3 quantization condition applied to the Hamiltonian Eq. (4.10), 
the zero-momentum states for each parity and their associated anomalous dimensions are 
given by: 



p 


states 


7 


-1 


W> & L/3 & L/ 3 ^lo) 


(5± s/7)\ 


1 


|0>, C/s^/slO) 


(2 ± V / 13)A 



It should be clear that the 1-loop spectrum of dimensions for closed operators with 
arbitrary values of L can be similarly obtained. 



4.1.2 Asymptotics 

Armed with the explicit expression Eq. (4.10), we can make some interesting observations 
about the pattern of operator dimensions for large classical dimension A^ ^ = 2L. In the 
limit L — > oo, the distribution of momenta becomes continuous, and the highest and lowest 
anomalous dimensions r y± for a given L can be approximated as elliptic integrals 

XL r 2w 

7± = XL ± — / dk^l0 + 6cos{k) 
47r Jo 

= (l± 4E ^ A ^ XL « (1 ± 1.542) XL. (4.16) 

Note that | '~y=i= I becomes arbitrarily large as L grows, so that dimensions of operators with 
different L begin to cross. Figure 2 shows the dimensions of generalized single trace operators 
of the form Tr(XYX . . . ) for the first few L and small values of A ~ 3N C /Nf — 1, and we 
can in fact see level crossing occur when A > 1/(1. 542L). 



4.2 Open Chains 

The case of open boundary conditions is slightly more complicated than that of a circular 
chain, but still exactly solvable at finite length. The solution was first written down by 
Dougot et. al. in 2004 [78], and clarified somewhat in [79]. We will study the Hamiltonian 

L L 

n=0 n=l 
L L 

= ^(4 + c„)(cJ l+1 - Cn+1 )-/iE(2ct Cn -l), (4.17) 

n=0 n=l 

acting on an open chain of length L + 2. Once again, we can diagonalize H with an 
appropriate Fourier transform and Bogoliubov transformation. However, we must take care 
with boundary conditions, since the spins at either end of the chain are nondynamical. In 
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Figure 2: Dimensions of flavor-singlet operators of the form Tr(XYX . . . ) in J\f = 1 SQCD in the 
Veneziano limit, as a function of A ~ 3N C /Nf — 1, up to 1-loop. 



particular, we should look for creation and annihilation operators that never mix eigenstates 
of 0q and 0f, +1 . Thus, let us consider the ansatz 

b[ = 4-dl k (4.18) 

L+X 

4 = Y.^ n ^n-Cn)+f{ky hn {c{ + C n ), (4.19) 
n=0 

where f(k) is to be determined, subject to the condition 

f(k)e ik{L+l) - f(-k)e~ ik{L+1) = 0. (4.20) 

Eq. (4.18) ensures that {b\, erg} = —i{b\,c\ — c } = 0, while Eq. (4.20) additionally ensures 
that = 0. In other words, b\ flips 0g and preserves c£ +1 , never mixing different 

boundary spins. Plugging our ansatz into the eigenvalue equation [H, b' k ] = e{k)b\, we find 

-2f{k){e- ik + h) = e{k) (4.21) 
-2(e ik + h) = f(k)e(k), (4.22) 

which has solution e{k) = 2^Jh? + 2hcos(k) + 1, as before. The boundary condition Eq. (4.20) 
becomes a quantization condition for the quasimomenta k, 

3i f< L + 2 ;> = -h. (4.23) 
sin(A;(L + 1)) v ! 
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When h > as is always case for us, Eq. (4.23) has L real solutions k E [0, 7r], 

along with an additional "bound state" solution of the form k = it + ij. Together, the b\ 
generate a Fock-space with dimension 2 L+1 , which exactly matches the number of states 
with fixed rightmost spin, say cr£ +1 = 1. Thus, after suitably normalizing our creation 
and annihilation operators, the Hamiltonian on this subspace can be written in the form 
Eq. (4.7). The Hamiltonian on the other subspace u z L+1 = — 1 is conjugate by an overall 
parity transformation and has an identical spectrum. 

Finally, let us identify which class of operators correspond to the highest and lowest 
energy states. When h *C 1, it's clear that the highest energy state has ferromagnetic 
boundary conditions o"q = <7f, +1 . One can show that no level-crossings occur as h varies, so 
that this is the case for arbitrary h > 0. The ground state is obtained by acting with all 
L + 1 lowering operators, and thus satisfies erg |0) = (— l) L+1 o-£ +1 |0). 



4.2.1 Example: Open Four-Field Operators 



Specializing now to SQCD, the open-chain Hamiltonian Eq. (3.11) can be written 

Kpen = A (l - ^ + A (blh ~ ^ VlO + ecos^), (4.24) 

with k subject to Eq. (4.23). As a simple example, let us consider open-chain operators with 
L = 1. Our quantization condition has two relevant solutions: 



_3 I y/13 \ / _3 _ a/13 
fc x = cos" 1 ( ^ — » 1.419, k 2 = cos' 1 ^ — ] « 7T - 1.087i. (4.25) 

Without loss of generality, we can consider states with positive rightmost spin. The 
highest and lowest energy states both have ferromagnetic boundary conditions, so they will 
be linear combinations of Q^XQ and Q^YQ. These correspond to the states 

|0> and 4X10), (4.26) 

which have the anomalous dimensions 



Q ± l - v/10 + 6 cos(A; 1 ) ± i ^10 + 6 cos(A; 2 )^ A = ^(1 ± v / 13)A. 



(4.27) 



Linear combinations of QXQ and QYQ correspond to the remaining states 

bl\0) and bl\0), (4.28) 
which then give the anomalous dimensions 

Q ± 100 + 6 cos(fci) T \ >/l0 + 6 cos(A; 2 )^ A = Q ± A. (4.29) 

Dimensions of operators with open flavor indices at larger values of L can be easily obtained 
in a similar fashion. 
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5 Conclusions 



In the present work we have shown that the 1-loop dilatation operator of M = 1 SQCD 
(in the electric Banks-Zaks limit) is equivalent to the 1-dimensional Ising spin chain in a 
transverse magnetic field, when acting on operators built out of scalar quarks. There are a 
number of directions that merit further study. One obvious direction would be to study the 
1-loop dilatation operator in other sectors of the theory, as well as higher-order corrections. 
In particular, one might ask whether the 2-loop dilatation operator in the scalar sector is 
solvable in any sense. It would be interesting to see if, for instance, a nontrivial ^-matrix 
for the Ising model fermions is induced at 2-loops. 

A related direction would be to study the dilatation operator in the magnetic dual 
description of M = 1 SQCD, which can be computed in perturbation theory when Nf ~ \N C . 
Because the dual description contains both dual quarks q a \ q a j and mesons M\, generalized 
single-trace operators containing scalars may be built out of chains of the color-adjoint 
(flavor-singlet) objects q(MM^) k q\ q^(M^M) k q, q(MM^) k Mq } and g t (M t M) fc M t g t . Since 
there is an infinite tower of possible states at each site, it would be interesting to understand 
if the 1-loop dilatation operator in this basis can be described by a simple noncompact spin 
chain, as well as to investigate its integrability properties. 

The transverse Ising model that we have encountered is one of the simplest examples of 
a spin chain with a quantum phase transition [75]. The transition from a paramagnetic to a 
long-range ordered phase occurs at h — 1, or naively when Nf/N c = 1, outside the conformal 
window and well outside the regime where perturbation theory is valid. Nevertheless, it is 
tempting to speculate that this phase transition persists in the full theory, but gets lifted to 
the bottom of the conformal window Nf/N c = 3/2. Studying higher-order corrections and 
also the dual magnetic description should help to shed light on whether or not this is the 
correct picture. 

Finally, there are many other 4D conformal field theories containing fundamental flavors 
that possess a Veneziano limit. Studying these theories is important because they help us to 
break outside of the M = 4 SYM universality class. One particularly simple example is Af = 2 
SQCD with Nf = 2N C flavors, where the 1-loop dilatation operator acting on scalars shares 
many features with that of M = 1 SQCD [68, 69]. Another example is non-supersymmetric 
QCD in the conformal window [39] and its many variations. It is an important goal of future 
research to better understand the extent to which these theories have solvable planar limits, 
and in particular to understand the role that supersymmetry plays in creating the kinds of 
simplifying features that have enabled so much progress in unraveling the structure of M = 4 
SYM. 
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